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Abstract. We analyze the static potential of a quark-antiquark pair at T > T c , where T c is a temperature of a deconfinement 
phase transition in QCD. We discuss the possibility that the non-perturbative part of this potential can be studied through the 
modification of the correlation functions, which define the quadratic field correlators of the nonperturbative vaccuum fields. 
We use the non-perturbative quark-antiquark potential derived in this way and the screened one-gluon-exchange potential 
with T-dependent Debye screening mass to calculate J/y, T and Sl bhh binding energies and melting temperatures in the 
deconfined phase of the full 2-flavors QCD. 
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INTRODUCTION 

There is a significant change of view on physical prop- 
erties and underlying dynamics of quark-gluon plasma 
(QGP), produced at RHIC, see e.g. [1] and references 
there in. On quite general grounds it is expected, that fun- 
damental forces between quarks and gluons get modified 
at finite temperature. Instead of behaving like a gas of 
free quasiparticles - quarks and gluons, the matter cre- 
ated in RHIC interacts much more strongly than origi- 
nally expected. Recall that above the deconfinement tem- 
perature T c the perturbative one-gluon-exchange poten- 
tial is expected to be exponentially screened at large dis- 
tances (r 3> l/T) [2]. Moreover, heavy quark free en- 
ergy F@@(r, T) shows large asymptotic value, F@Q(°°,T) 
[3]. This value can be explained only by non-perturbative 
(NP) effects, since perturbative one-gluon-exchange po- 
tential, even with increased <X s (r), cannot produce sim- 
ilar effect. Therefore it is more appropriate to describe 
the (NP) properties of the QCD phase above T c in terms 
of the NP part of the QCD force rather than a strongly 
coupled Coulomb force. 

At T = 0, the non-perturbative quark-antiquark poten- 
tial is V np (r) = or, where a is the 51/(3) string ten- 
tion. This potential has been used extensively in poten- 
tial models. At T > T c , a = 0, but that does not mean 
that the NP potential disappears. Attempting to guess the 
form of the non-perturbative potential we address to the 
Field Correlator Method (FCM) (see [4] and references 
there in). Within FCM the NP QQ potential above T c 
was suggested to occur due to the non zero correlation 
function D\ (x) that is one of the two functions which de- 
fine the quadratic field correlators of the nonperturbative 
vaccuum fields . The most direct prediction of this ap- 
proach is the existence of bound states of heavy cc and 



bb mesons and heavy bbb baryons above T c . In a recent 
paper [5] we calculated binding energies for the lowest 
QQ and QQQ eigenstates (Q — c,b) using both the NP 
potential of the FCM and the screened Coulomb poten- 
tial with the temperature dependent Debay radius calcu- 
lated in pure 51/(3) glue theory. In this talk we discuss 
and extend the results of this analysis, in particular, by in- 
cluding the effect of Debye screening in the full 2-flavor 
QCD. 



FCM AT FINITE TEMPERATURES 

The approach is based on the study of the quadratic field 
correlators < trF )iV (x)^(x, 0)Fx a (0) > (x is Euclidian), 
where 4>(jc,0) is the Schwinger parallel transporter nec- 
essary to maintain gauge invariance. These correlators 
are expressed in terms of two scalar functions, D(x) and 
D\{x) [4]. At T = 0, the string tension a is expressed 
only in terms of D(x): 



a = 2 J dX J dvD(y / X 2 + v 2 ). (1) 



At T > T c one should distinguish between electric and 
magnetic correlators D E (x), D H (x), Df(x), and Dj (x), 
and, correspondingly, between a E and a . It was argued 
in [6] and later confirmed on the lattice [7] that above the 
deconfinement region D E (x) and a E vanish, while the 
colorelectric correlator Df (x) and colormagnetic corre- 
lators D H (x) and D f (x) should stay unchanged at least 
up to T ~ 2 T c . The correlators D H {x) and Df (x) do not 
produce static quark-antiquark potentials, they only de- 
fine the spatial string tension (7 V = a H and the Debye 



mass mo « y/o^ that grows with the temperature in the 
dimensionally reduced limit [8]. 

The NP static QQ potential at T > T c originates from 
the color-electric correlator function Df (x) 

l/T r 

V np (r,r) = j dv{\-vT)j XdXD*{x). (2) 



In the confinement region the function Df (x) was calcu- 
lated in [9] exploiting the connection of field correlators 
to the gluelump Green's function 1 

Dfw=B !^p(dM ! (3) 

where B = 6a/oyM , a/ being the freezing value of the 
strong coupling constant, oy is the sting tension at T = 0, 
and the parameter Mq has the meaning of the gluelump 
mass. Above T c the analytical form of D\ should stay 
unchanged at least up to T ~ 2T C . The only change is 
B — > B(T) = % (T)B, where the T-dependent factor 

§(r)=(l- 0.36^1^) (4) 

is determined by lattice data, see Eq. 52 of Ref. [11]. 
Substituting (3) into (2) and integrating over X one ob- 
tains V„ p (r,T) = V(oo,T) - V(r,T) where 



V{eo,T) = 

and 

V(r,T)-- 




o 



One observes the characteristic feature of the static po- 
tential V np (r,T) produced by the correlator Df (x): the 
potential gives rise to the constant term V(°°,T) in the 
QQ interaction at large distances, which can be viewed 
upon as the sum of selfenergies of Q and Q. 

We can now exploit the relativistic Hamiltonian tech- 
nics [13] successfully applied for mesons, baryons, glue- 
balls and hybrids in the confinement phase. This tech- 
nic does not take into account chiral degrees of free- 
dom and is applicable when spin-dependent interaction 
can be treated as perturbation. Therefore below we con- 
sider only heavy quarkonia and heavy baryons, leaving 



TABLE 1. Parameters of the quark- 
antiquark potentials in units of GeV. 
Both sets of parameters correspond to 
K(oo,r c )= 0.505 GeV 



tlf 


T c 


M 


m d (T c ) 





0.275 


0.9 


0.793 


2 


0.165 


1.08 


0.545 



light quarkonia with chiral symmetry restoration to an- 
other publication. 

Recall that, in the framework of the FCM, the masses 
of heavy quarkonia are defined as 

2 

M QQ = 77^+^6+ E (™q,Hq), (7) 

Eo(mQ,jj,g) is an eigenvalue of the Hamiltonian 

H = H Q + V np + V 0G E, (8) 

where we have omitted spin-dependent and self-energy 
terms proportional to 1/ [Iq. In Eq. (8) Voce is the one- 
gluon-exchange potential which is expected to be expo- 
nentially screened at large distances 

VoGE(r,T) = -^exp(-m d (7», (9) 

m d {T) being the Debye mass. In Eq. (7) mq are the bare 
quark masses, and einbeins jj,g are treated as c-number 
variational parameters 2 . With such simplifying assump- 
tions the spinless Hamiltonian takes an apparently non- 
relativistic form, with einbein fields playing the role of 
the constituent masses of the quarks. In what follow we 
take m c —1.4 GeV, = 4.8 GeV. As in the confinement 
region, the constituent masses \Iq only slightly exceed 
the bare quark masses mg that reflect smallness of the 
kinetic energies of heavy quarks. The dissociation points 
are defined as those temperature values for which the en- 
ergy gap between V(°°, T) and Eq disappears. 

QUARK- ANTIQUARK STATES 

The non-perturbative quark-antiquark potential is de- 
fined by the two parameters B and M . In what fol- 
lows we take Of = 0.18 GeV 2 and a/ = 0.6, so that 
B = 0.648Mo, and vary the gluelump mass Mq around 
the central value Mq = 1 GeV in order to maintain 
the asymptotic value V(°°,T C ) = 0.505 GeV. This value 



1 Recall that gluelumps are actually bound states of the gluon field in 
a static color-octet source that have been studied first in Lattice QCD 
[10]. 



2 The eigenvalues Eo(mQ, [Iq) are found as functions of the bare quark 
masses mg and einbeins Hq, and are finally minimized with respect to 
the )1q. Once mq is fixed, the quarkonia spectrum is described. 



TABLE 2. J/y states above the deconfinement region. All the quantities except 
for T/T c and ro are given in units of GeV, the dimension of ro is GeV -1 . 



T/T c 


m d 


v(«o,r) 




Eo(T)-V(«o,T) 






1 


0.545 


0.505 


1.462 


-0.026 


6.89 


3.281 


1.2 


0.609 


0.433 


1.438 


-0.0098 


11.45 


3.334 


1.3 


0.640 


0.398 


1.426 


-0.0046 


16.86 


3.164 


1.4 


0.671 


0.365 


1.415 


-0.0013 


25.51 


3.164 



TABLE 3. T states above the deconfinement region. The notations are the same 
as in Table2. 



T/T c 




V(oo,r) 


Mi 


£b(r)-v(oo,r) 


ro 


M bb 


1 


0.545 


0.504 


4.948 


-0.345 


1.17 


9.768 


1.6 


0.733 


0.302 


4.954 


-0.182 


1.54 


9.725 


2.0 


0.853 


0.187 


4.937 


-0.102 


2.01 


9.688 


2.8 


1.082 





4.837 


-0.008 


7.88 


9.592 



TABLE 4. Dissociation temperatures 
(in units of T c ) for cc, bb, and Slbbb 
states. Q. ccc is unbound both for nj = 
and «/ = 2. 



IS bottomoniium undergos very little modification till 
T ~ 2T C . The melting temperatures for the J/y/ and T 
are shown in Table 4. 





J/Y 


T 


&bbb 




n f = 
n f = 2 


1.29 
1.48 


2.57 
2.96 


1.8 
2.35 


QQQ BARYONS 



agrees with lattice estimate for the free quark-antiquark 
energy 3 . The strong coupling constant was taken a s = 

0. 35 4 . The parameters of the potential are listed in Table 

1, where for the reference we also indicate the values of 
the Debye mass m^(r c ) [12]. 

Some details of our calculation for the full n/ = 2 
QCD can be inferred from Tables 2, 3 5 . At T = T c we 
obtain the weakly bound cc state. The melting tempera- 
ture is <~ l.3T c for rif = and 1.48 T c for rif = 2, see 
Table 4. The charmonium masses lie in the interval 3.1 - 
3.3 GeV. Note that at the melting point r{J/y) oo that 
is consistent with nearly-free dynamics. 

As expected, the T state is much more bound and 
remains intact up to the larger temperatures, T <~ 2.3 T c . 
This is in agreement with the lattice study of Ref. [15]. 
The masses of the L = bottomonium lie in the interval 
9.6-9.8 GeV, about 0.2-0.3 GeV higher fhan_9.460 GeV, 
the mass of T(15) at T = 0. At T = T c the bb separation 
ro is 0.25 fm that is compatible with ro = 0.28 fm at 
T = (at the melting point ro — > °°). Note that the 



3 However, the difference in the parameter Mq causes the small differ- 
ence of V(t», T) for T > T c . 

4 The account of the running Ct s (r) in the Coulomb potential produces 
a tiny effect as compared with the case of a constant a s = 0.35 both for 
the energies and wave functions [14]. 

5 The corresponding results for the pure gluodynamics (rif =0) are 
given in Ref. [14] 



The three quark potential is given by Vqqq = 
2 £i<j ^QQ{ r ij^T), where j is the color factor and Vqq 
is the sum of the perturbative and NP quark-antiquark 
potential. We solve the three quark Schrodinger equation 
by the hyperspherical harmonics method. The wave 
function in the hypercentral approximation is written as 

1 u(R,T) 



*V(R,T) 
where the hyperradius 



3 /?5/2 



(10) 



(11) 



is invariant under quark permutations. Averaging the 
three-quark potential over the six-dimensional sphere 
one obtains the one-dimensional Schrodinger equation 
for the reduced function u(R,T) 



d 2 u{RJ) 

dR 2 
where V(R,T 



Eq 



15 

oge(R, T) 

n/2 



-f(R,T) 



u(R,T)=0, (12) 



y„ p (R,T) and 



y OG E(R,T) = - 



16 a s 
1% 



exp(-m d (T)R) . 2 



R 



sin 2 (20), 
(13) 



y np (R,T)=V(oo,T)-W^-x 



7ZMq 



k/2 

I 



RKx(R) 



M 



e- R (l+R) )sm 2 (2G)d9, (14) 



TABLE 5. Clbbb state above the deconfinement region. The interquark distances 



T/T c 




r(oo,r) 


Mi 


£b(r)-r(»,r) 


V<R 2 > 


Mbbb 


1 


0.545 


0.757 


4.962 


-0.327 


3.44 


14.837 


1.4 


0.672 


0.548 


4.926 


-0.185 


4.22 


14.768 


2.0 


0.853 


0.281 


4.919 


-0.192 


7.56 


14.641 


2.3 


0.940 


0.166 


4.830 


-0.0034 


18.50 


14.563 


2.4 


0.969 


0.131 


4.812 


+ 0.0021 


32.40 


14.533 



V (oo, T) being given by Eq. (5), and R = 2M Rsin 9 /jUg. 
In Eq. (14) we use the approximate expression for the 
non-perturbative QQ potential (6) 

y(r ' r) * { xKl(x) ~ Wo expi ~ x){1 +x) ) ' 

(15) 

where x = M§r and K\ (x) is the McDonnald function, 

The temperature dependent mass of the colorless QQQ 
states is 

3 m 2 n 3 

M QQQ = nlT + o^e+ £ o(m e ,Al e ), (16) 

where jiig are now defined from the extremum condition 
imposed on Mqqq in (16). The bound QQQ state exists 
ifE (m Q ,H Q ) < Vqqq(<*>,T), where 

fQQQ^J) = \v(oo,T). (17) 

There is no bound £2 CCC states 6 but the survives up 
to T ~ 1.8 — 2.4 7^ (depending on n f), see Tables 4, 5. 

CONCLUSIONS 

The static QQ potential has been extensively investi- 
gated within the FCM and provides useful tool to study 
in-medium modification of inter-quark forces. This po- 
tential provides also useful quantitative insights into the 
problem of quarkonium binding in QGP. In particular, 
the color electric forces due to the nonconfining corre- 
lator D\ survive in the deconfined phase and can sup- 
port bound states at T > T c . In this paper, we used a 
FCM approach to the problem of heavy quark potentials 
at finite temperature. We have calculated binding ener- 
gies and melting temperatures for the lowest eigenstates 
in the cc, bb, and bbb systems neglecting spin-dependent 
and self-energy terms in the Hamiltonian. We find that 



the ground state of J/y survives up to T <~ 1 .3 — 1 .5 T c , 
and there is no bound Cl ccc state at T > T c . On the other 
hand, the bb and bbb states survive up to higher temper- 
ature, T ~ 2.6 - 3.0 T c and T ~ 1 .8 - 2.4 T c for n f = 0, 2, 
respectively. The results suggest that the systems are 
strongly interacting above T c . 
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